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Quantum-mechanical calculations of muon transfer between muonic hydrogen and an oxygen 
nuclei for s waves and collision energies in the range 10 -3 - 10 3 eV, are presented. Close-coupling 
time-independent Schrodinger equations, written in terms of hyperspherical elliptic coordinates were 
integrated along the hyper-radius to obtain the partial and total muon-transfer probabilities. The 
results show the expected Wigner-Bethe threshold behavior up to collision energies of the order of 
10 -2 eV and pronounced maxima at 10 2 eV which can be interpreted in terms of crossings between 
potential energy curves corresponding to the entrance channel state (fJ.p)is + O and two product 
channels which asymptotically correlate to p + (0^) n =5,6- The population of the final states with 
different orbital angular momenta is found to be essentially independent of energy in the range 
considered in this work. This can be attributed to a strong selection rule for the conservation of the 
quantum number associated to one of the elliptic hyperangles. 

PACS numbers: 36.10.Dr 
Keywords: 



I. INTRODUCTION 

Negative muon transfer between muonic atoms 
(muonic hydrogen, for instance) and other atoms or 
molecules has been extensively studied in the framework 
of muon catalyzed nuclear fusion (see Ref. [l[ and lit- 
erature cited therein). Also, the structural and spec- 
troscopic properties of these species are of interest for 
metrology and in tests of quantum electrodynamics @, Q . 

Recently, several theoretical @, 0, H, H, 0] and ex- 
perimental Q works have considered the problem of 
muon transfer from the muonic hydrogen to an oxygen 
molecule. Since the muonic hydrogen has to approach 
one of the oxygen nuclei very close in order for the muon 
to be transferred [9j], the process can be described as 
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Although there have been several full three- 
dimensional calculations of muon transfer rates at 
low energies between muonic-hydrogen and \ow-Z atoms 
(see literature cited in Ref. Q), there is none when 
the transfer involves nuclei with Z > 3 . Indeed, as Z 
increases there is a larger initial-channel polarization 
and a stronger final-channel Coulomb interaction which 
make the full quantum calculation computationally 
heavy. Thus up to now only approximate calculations 
have been performed for the muon-transfer rate between 
muonic-hydrogen and oxygen [J, @, 0, [IfJ- We present 
here the first numerically converged three-dimensional 
calculations for reaction ([1]) for s waves and collision 
energies in the range 1CP 3 — 10 3 eV. Since for the 
entrance channel the centrifugal barrier for J = 1 is 
about 0.1 eV, the calculations presented here can be 
considered as full 3-dimensional up to thermal energies. 



The calculation were performed as follows. Hyper- 
spherical elliptic coordinates [III EH have been used. A 
piecewise diabatic basis set on the hyperspherical angles 
was used to expand the wave function. The resulting 
close-coupling time-independent Schrodinger equations 
in the hyper-radius were solved using a de Vogelaere al- 
gorithm and the partial and total muon-transfer proba- 
bilities were determined by the standard S-matrix anal- 
ysis at large distances. Since for energies below 10 _1 eV, 
the muon-transfer process studied here is equivalent to 
an ultra-cold collision (de Broglie wavelength, A > 1 A, 
much larger than the effective range, a ~ 0.1 A, of the 
potential interaction) , special care had to be taken to the 
asymptotic analysis in the entrance channel. 

The paper is organized as follows. Section II intro- 
duces the model and the methodology used in the calcu- 
lations. Section III presents the calculated muon-transfer 
probabilities together with their interpretation in terms 
of simple Landau- Zener and threshold models. Finally, 
section IV is devoted to the conclusions. 



II. METHODOLOGY 

We start with the two mass-scaled Jacobi sets of co- 
ordinates (.Ri,T"i) and (i?2,'"2) adapted to the entrance 
and exit channels of reaction |T]), respectively (see fig. []}. 
They are defined by 
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In spherical coordinates, the system can be described 
by three Euler angles specifying the overall orientation, 
the two distances (iij,r,*) and the angle 7, between the 
two vectors Ri and r i (see figure [1} . Two sets of Delves 
hyperspherical coordinates are then defined by the com- 
mon hypcr-radius 



the hyper-angles 



(8) 



FIG. 1: The two sets of mass-scaled Jacobi coordinates cor- 
responding to the entrance and exit channel of reaction |T}. 
The center of mass positions are not in scale. 



where the set (i = l,j = 2, k = 3) corresponds to (O, p, 
p), and (i = 2,j = 3, k = 1) to (p, p, O). The reduced 
masses rtiijk and m,j t k in Eqs. ^ and Q, are given by 



rrii {rrij + m k ) 
rrii + rrij + m^ 

while m is choosen as 



mi rrij rrik 



rrij rrik 
rrij + m k 



1/2 



, (4) 
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v rrii + rrij + rrik , 
These sets are related by the orthogonal transformation 
— cos 9,, — sic 
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where 



tan 6^ = 



m 



(6) 
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giving in our case 9 U ~ 19° 



tanfo/2) = -f ; i = l,2. 
Hi 



(9) 



and the 7, angles. 

The relationship between the two sets of Delves shape 
angles is given by 

COSX2 
sinx2 cos 72 
sinx2 sin 72 

cos(20 M ) sin(2^) 0\ / cosxi 
-sin(2^) cos(20 M ) sin X i cos 7l |(10) 
1 J \ sinxi sin 71 

The hyperspherical elliptic coordinates are defined by 

V = Xi ~X2, -20p<ri<20„ (11a) 
£ = X1+X2, 2^<£<2tt-2^, (lib) 

with the volume element 

_2 



dr = p dp 



7T 



4sin(20 M ) 



(cos(»7) - cos(£)) dr] d£. (12) 



In terms of these coordinates the kinetic energy oper- 
ator for total angular momentum zero, is given by 



T = — 



1 d 



d 16 



2m \ p 5 dp^ dp p 2 cos(rj) — cos(£) 



- ^cosfa) - COB(2^)J - - - (co^) - cos(20 M )J - 



(13) 
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For a given value of the hyper-radius /?, the total wave 
function ?/>(p, ??, £) is expanded in terms of a basis set of 
N c h functions 4>i(r], £; p) depending on the hyperspherical 
angles 77 and £. We use a diabatic- by-sector representa- 
tion. In each sector p n — Sp n < p < p n + Sp n ;n = l,...,N p 
we write: 



where ^(77, £; p n ) are eigenstates of the Hamiltonian at 
fixed p n distances. Their calculation requires the solution 
of a bound state problem for the Coulomb potential 



V 



8e 2 



tCo X, 



\X U Xr 



\X V Xr 



(15) 



4>(PiV,0 = -573 F i(fi)<l>i(V>tlPn) 



presenting two attractive singularities at (7/, £) = 
(74) (±29^,29^) corresponding to a vanishing muon-oxygen 
and muon-proton distance. 
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This bound state problem can be rewritten as: 

f i6h 2 - - 1 

(16) 

where 

d d 
L(u) = — (cosu-cos(2^)) — (17) 

and 

W(M,0 = [COS(»J)-C08(0] [^(Pn,»7,0-Ci(Pn)]- ( 18 ) 

Equation (|16p can be viewed as a zero eigenvalue prob- 
lem depending parametrically on the potential [cos (77) — 
cos(£)] [V(p n ,r], £) — Ci(p n )], in the sense that the constant 



£i(pn) is being adjust in such a way that the operator on 
the left hand side of equation (flU)) has zero eigenvalue. 
The renormalized potential W(fl n ,T], £) has two impor- 
tant features. One is that the two Coulomb singularities 
are regularized, the other being that it is approximately 
separable. We can therefore write 

W(p n , 77, = W v (p n , rj) + W^(p n , + AW( Pn , V , 0, 

(19) 

where W n (p n ,rj) = W(p n , 77, 2 M ) and W^(p n ,0 = 
W(p n , -26^,0 -W( Pn , -20^,29^). 

The 2-dimensional problem can be solved by exploit- 
ing this approximate separability of the potential. For 
instance, defining 77 = 2 9^f) with — 1 < fj < 1, we get 
from equation (fl6|) 
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W v {p n ,rj) 



Vk{mPn) = t'kHpn) Vk{m Pn) 



(20) 



The similarity between the differential operator in ([20]) 
(in particular in the limit 9^ — > 0) with the one defining 
Legendre polynomials, suggests to use the latter as basis 
set functions for expansion of tpkiVi p n ). We solve a simi- 
lar problem for £ using the W^(p„, £) potential and obtain 

Pn) eigenfunctions and the e^(/9„) eigenvalues. We 
then iterate (using a bi-section method on e{) until we get 

4. + e \ (pn) = 0. Once the separable basis set is 
obtained, solutions of (fl6|) for the full non-separable po- 
tential are obtained by diagonalizing the representation 
matrix of the full Hamiltonian in the product basis. In 
Fig. [2] we present the calculated energies ej as a function 
of the hyper-radius p. The origin of energies has been 
chosen to coincide with the asymptotic limit of the en- 
trance channel (p p) n =i + O. The calculations presented 
in this work cover the energy range between this limit 
and the p + (p O) n= io threshold at about 1 keV. 

The N c h coupled equations are integrated along the 
hyper-radius p using the de Vogelaere algorithm [l3j |. 
This provides a logarithmic derivative matrix Z at pm = 
Pn p + &n„- For the energy range considered here, we in- 
cluded 88 channels: ((p/u) n =i-2 + andp+(/jO) n= i_g). 
The integration of the coupled equations was performed 
from the origin to p en( ^ ~ 200 ~ 1 A. The asymp- 
totic analysis has been performed using the appropriate 
Jacobi coordinates for the entrance and for the product 
channels. Elementary asymptotic wave functions for the 
different final arrangement channels A = 1, 2 are written 
as products of translational functions f n i(Rx), Coulomb 
bound wave functions C n i(r\) and normalized Legendre 
polynomials: 

Vwa(-Ra,t-a,7a) = fne(Rx) C n i{r\) P t (cos 7*). (21) 
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FIG. 2: Calculated e, energies of the hyperspherical elliptic 
basis as a function of the hyper-radius. 

We use these elementary asymptotic wave functions to 
form N c h physical solutions whose forms are given by 

+ ^rcW(#V,rv,7A')^'A'^«n (22) 

n'l'X' 

where the superscripts S and C refer to the type of 
translational functions (S for sine- type functions for open 
channels and exponcntiallly decaying functions for closed 
channels, C for cosine-type and exponentially growing 
functions, respectively). If a Coulombic interaction is 
still important (as for instance in the product channels 
Op+p in our problem), Coulomb wave functions are used 



4 



instead. 

Wc define the matrices F and F' as the projections 
of the elementary asymptotic functions on the fa hyper- 



J 



spherical basis. This projection is performed at the max- 
imum hyperspherical radius pm 



Fi,n£\ = (fa\lpnl\) p=p„ 



dr)d£, [cos(?7) - cos(£)] fa{v, & Pm) ip n ex(Rx, r\, J\)\ f 



F' 



■ ex 



\dip n ex/dp)p=p M = / dr]d4 [cos(rj) - cos(£)] fa(v,^,PM) dip n e\(Rx > r\,j\)/dp\ 



PM 



(23a) 
(23b) 



The N c h linearly independent solutions which result from 
the propagation steps up to pm are linear combinations 
of the N c h asymptotic solutions given by equations (f2"3"| . 
This can be restated as an equality of the logarithmic 
derivative matrix Z in the hyperspherical basis 



(F' S + F' C K)(F i 



F c K)~ 



(24) 



The K matrix, and then the S matrix, can be extracted 
from (EH. 



III. RESULTS OF THE CALCULATIONS 

We have performed calculations of the reaction 

(pn) n =i+0 8+ ^p+( P 0) 7 n + (25) 

for collision energies in the range 10 -3 — 10 3 eV and for 
total angular momentum J = (s waves). For large 
R distance between the oxygen nuclei and the proton, 
the potential for the entrance channel of reaction |25|) 
behaves as V = —aZ 2 e 2 /2R 4 , where Z = 8 and a = 
(9/2)(?i 2 /m M ,pe 2 ) 3 . Since with this potential the centrifu- 
gal barrier has a height of (h 2 J{ J +1) /m ^ pil Ze) 2 /Ha, for 
energies below 0.1 eV the partial wave J = is the only 
one which contributes to the cross sections. Thus for 
thermal energies these calculations are essentially exact 
full 3-dimensional. 

Fig. [3] presents the total probabilities for muon trans- 
fer as a function of the energy, as well as the partial 
muon-transfer probabilities into the p + (p 0) n= e and 
p+{p 0) n= 5 channels. The other open channels have ncg- 
ligeable populations. This is in complete agreement with 
what was found in approximate calculations @. These 
results can be qualitatively understood by inspection of 
figure [2j Starting in channel (pp) n= i + O, the system 
crosses diabatically the channel p + (pO) n= j. Muon 
transfer is completely negligible as the coupling is very 
small compared with the collision energy. The couplings 
to channels p + (p 0) n= 6 and p + (p 0) n= s are larger 
as evidenced by avoided crossings. The other channels 
p + (p 0)„<5 are weakly coupled to the initial one and 
they are not expected to be populated significantly. 





0.6 r 




0.5 




0.4 












0.3 












0.2 




0.1 










| total | 




















|» = 8| 









0.001 0.01 



0.1 1 

E (eV) 



10 



100 



1000 



FIG. 3: Partial and total muon-transfer probabilities for re- 
action (|25[) as a function of the collision energy. 



Actually, a simple calculation involving only three 
channels ((p/z)„=i+0, p+{pO) n =5-6,i=o) and using the 
semi-classical Landau-Zener formula for the two cross- 
ings provides qualitative agreement with the exact re- 
sults down to energies of the order of 1 eV (see figure [4]). 
For lower energies, quantum threshold behavior domi- 
nates the semi-classical Landau-Zener approximation and 
for very low energies Wigner's threshold laws predict a 
E 1 / 2 dependence. This low energy region corresponds 
to ka <C 1, with k = (2m ow £) 1 / 2 /'i and a being the 
range of the potential in the entrance channel. Defin- 
ing the effective range of the potential by V(a) < E, we 
get E <C 0.1 eV. This is exactly what it is found in our 
calculations (see figure El • 

In figures [5] and [HI we present relative populations of 
Op n= 5-6.g levels for collision energy E = 0.1 eV. This 
distribution is almost independent of collision energy in 
the energy range considered here. This fact and the par- 
ticular form of the distributions in figures [5] and [5] can 
be understood as follows. From an inspection of figure 
[21 it is clear that only one of the set of states correlating 
asymptotically to a given Op n manifold is significantly 
coupled to the initial pp n= \ state. Due to the approx- 
imate separability of the potential, each of these states 
can be assigned to approximate hyperspherical elliptic 
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FIG. 4: Total muon-transfer probabilities for reaction (|25fl as 
a function of the collision energy using a simple three channel 
problem and Landau-Zener semi-classical approximation. 
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FIG. 5: Exact (full line) and approximate (dotted) relative 
populations of the I levels for n = 5 and collision energy 
£ = 0.1eV. 



quantum numbers (n n , n^). Since the initial p/j, n= i chan- 
nel is colinearly dominated and > corresponds to 
excitation away from the colinear configuration, the fi- 
nal states significantly coupled to the initial channel sat- 
isfy the approximate selection rule: = 0. At large 
distances, hyperspherical coordinates tend gradually to- 
ward parabolic coordinates [l^- The final (n v ,n^ = 0) 
states can therefore be identified to (h v ,n^ = 0) states 
with = n — 1 and = n^). The final approximate I 
distribution is then given by the square of the projection 
of this states onto the final (n, t) spherical basis set with 
the result [3]: 

P a J cx \C((n - l)/2, (n - l)/2, -n/2, -n/2; £, 0)| 2 (26) 

where C stands for a Clebsh-Gordan coefficient. We have 
represented this approximate distribution in figures[5]and 
[6l Clearly, the — selection rule is a very good ap- 
proximation. 

Experimentally only the populations of the fine- 
structure levels j are eventually measurable. Thus, from 
the orbital angular momentum I probabilities one gets 

p n,j oc 21 + \ ^ \C{t,l/1>™>t,m,j - m e ;j,mj)\ 2 

£ rrij ,mf 

(27) 

The results are presented in figure [7] 

IV. CONCLUSIONS 

We have presented 3-dimensional calculations of muon 
transfer probabilities between muonic hydrogen and oxy- 
gen for relative translational energies between 10~ 3 and 
10 3 eV and total angular momentum J = (s waves). 
Since the centrifugal barrier for J — 1 is of the order of 
0.1 eV, for thermal and lower energies the calculations 
are essentially exact. 
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FIG. 6: Exact (full line) and approximate (dotted) relative 
populations of the £ levels for n = 6 and collision energy 
E = 0.leV. 
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FIG. 7: Relative populations of the j levels for n — 5 (dotted) 
and n — 6 (full) for a collision energy E = 0.1 eV. 
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The results for the total and partial probabilities can 
be interpreted qualitatively by a simple three channels 
model and Landau-Zener semi-classical approximation 
for energies between 1 and 1000 eV. For energies below 
0.1 eV, the probabilities follow the expected Wigner's 
threshold laws behavior. 

The fine-structure levels populations are essentially in- 
dependent of energy and almost independent of the prin- 
cipal quantum number n. This has been attributed to an 
approximate selection rule for the elliptic hyperspherical 
quantum number nj. 

In conclusion, the hyperspherical elliptic coordinates 



introduced by Tolstikhin et al |ll|, are particularly well 
suited to the treatment of this type of reaction. 
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